Quantum state swapping via qubit network with Hubbard interaction 
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We study the quantum state transfer (QST) in a class of qubit network with on-site interaction, 
which is described by the generalized Hubbard model with engineered couplings. It is proved that the 
system of two electrons with opposite spins in this quantum network of N sites can be rigorously 
reduced into N one dimensional engineered single Bloch electron models with central potential 
barrier. With this observation we find that such system can perform a perfect QST, the quantum 
swapping between two distant electrons with opposite spins. Numerical results show such QST and 
the resonant-tunnelling for the optimal on-site interaction strengths. 
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Introduction. For implementing quantum information 
processing based on the scalable systems, the solid-state 
data bus is a necessary element to coherently integrate 
two or more qubits and transfer the quantum informa- 
tion among them 0. Recently increasing investigations 
have explored the possibilities to transfer quantum states 
through a class of solid-state data buses, the artificial spin 
chain with engineered nearest neighbor (NN) couplings 
0. El Some novel physical mechanisms 
have been discovered behind the protocols of quantum 
state transfer (QST) based on the spin chain systems. 
For example, it is discovered that the gap structure of 
spectrum of the strongly correlated systems is responsi- 
ble for the role of data bus Q ; and the spectrum march- 
ing parity symmetry is a sufficient condition for perfect 
QST HQ. 

We also notice that most of the explorations for QST 
are carried out only for the non-interacting systems or 
single-particle quantum states. And it seems that the on- 
site Coulomb interactions may destroy the quantum co- 
herence of transferred state. In this letter, we will study 
the influences of on-site interactions on the dynamic pro- 
cess of QST by making use of the generalized Hubbard 
model with engineered NN couplings as same as that in 
the artificial Bloch electron model in ref. Q . 

To give prominence to our central context we only 
consider the simplest interacting system with only two 
electrons of opposite spin involved. The main result we 
achieved is the discovery of the novel model reduction 
that the iV-site two-electron engineered Hubbard model 
can be decomposed into N single-particle engineered 
models on Z-sites chain (I = 1, 3, 2N — 1), but with an 
additional central potential barrier (CPB). This discov- 
ery enlightens us to conjecture the possibility of imple- 
menting the perfect quantum information swapping since 
the reduced models still keep the mirror symmetry. The 
detailed numerical simulations demonstrate there indeed 
exists such perfect quantum state swapping even certain 
on-site repulsion U is considered. 

Engineered Hubbard model and its reduction. Our 
model for quantum state swapping is the generalizations 
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FIG. 1: (color on line) (a) The engineered Hubbard model of 
N sites with two electrons, (b) The two-leg ladder of spin- 
less Bloch electrons, which is equivalent to the above Hubbard 
model. The rung represents the on-site interactions. We label 
the sites in the two legs with a standard angular momentum 
basis \JM) in two opposite orders; (c) The above equivalent 
two-leg network can be further reduced according to the prod- 
uct representations SO(3) ® 5*0(3), as the direct sum of N 
central potential barrier model. 



of the engineered spin model in ref. |j| (see the Fig. la) 
by adding the on-site Coulomb interactions. It can also 
be regarded as an engineered Hubbard model [lOj with 
artificial hopping. The model Hamiltonian reads 
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where cj g is the creation operator of electron at site j 
with spin a =f , J, and U is the on-site repulsion. The 
hopping integral is engineered as tj — sj j(N — j). It 
has been widely studied in connection with correlation 
effects in narrow-band solids and the concept of entan- 
glement ^H- If the on-site interaction is absent, i.e., 
U = 0, it has been shown in 0. El that an arbitrary 
many-particle state can be transferred to its mirror coun- 
terpart perfectly after the time r = ir/2. It is due to the 
fact that the energy-level structure and the parity of the 
corresponding cigenstate satisfy the spectrum-symmetry 
matching condition (SSMC) introduced in |9(: Let cf> n be 
the common eigen-function of H and mirror operator R 
with the eigen- values e„ and p n respectively. It is easy to 
find that any state ip at time r can evolve into its symmet- 
rical counterpart Rip if the eigenvalues e n and p n match 
each other and satisfy the SSMC exp(— ze„r) = p n . 
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It is easy to imagine that the energy levels should be 
shifted by ceratin deviations from the original spectra 
when the on-site repulsion is switched on. Nevertheless 
there still exists the possibility that the new set of shifted 
energy levels satisfy the SSMC for an appropriate U since 
the certain symmetry remains as will be shown in the fol- 
lowing discussions. We first illustrate our analysis along 
this direction schematically in the Suppose that there 
are only two electrons with opposite spin in the engi- 
neered Hubbard model (see Fig. la). The on-site inter- 
action occurs only when the two electrons occupy a same 
site. Alternatively, the Hubbard chain is equivalent to 
the spinless Bloch electron network with two legs and N 
rungs (see Fig. lb). The on-site interaction is denoted 
by the rungs. Each site in the leg corresponds to a sin- 
gle electron Bloch state, \j)^ — |0) or \j)^ = |0), 

./ i. •->.••• ..v. 

According to the ref. ji| , we can associate these states 
to the angular momentum states 

|J,M) t = \ J + M + l) v \J,M) l = [J-M + l)^ (2) 

where a given J = (N - l)/2, M = J, J - 1, • • • , - J + 
1, — J. Then it is easy to check that for the engineered 
couplings, the lowing operator of angular momentum can 
be realized in terms of the fermion operators as 

& - 4 t) -^ t) = E^4t c ™ 0) 
j 

which generates the group 50(3) together with = 

(J^V and ji a) = Yt,j<T3 c j,a c j,(T where a =|, \. Then we 
can rewrite the Hamiltonian as 

H = 2<4 T) + 2,4 X) + V; 

V = U^2\J,M;J,-M)(J,M;J,-M\, (5) 

M 

where the two-particle associated state \J,M; J, —M) is 
defined by | J, M; J, M') = \J,M)^®\J, M') ; . 

The intrinsic dynamic symmetry of the above gener- 
alized model is described as SO(3) ® SO(3) . Thus the 
addition theorem for two angular momenta [l2T | can be 
employed to reduce the representation of this general- 
ized Hubbard model according to the decomposition of 
the product representation ® = Y?l=o® d[L \ 
where is an irreducible representation of SO(3). The 
key point in our treatment is to expressed the on-site 
interaction term V as the sum of irreducible tensor op- 
erators. To this end, we use the Clebsch-Gordan coeffi- 
cients C^^.j^ = ((JJ)L,JM | J, Mi; J,M 2 ) to write 
the eigenvector of the total angular momentum \L, M) = 
\(JJ)L,JM) = Z Ml+ita=M Cjf$ 1 .j iMa \J,M 1 ;J,M 2 ). 



From the corresponding inverse transformation, the in- 
teraction term can be decomposed as 

V = UJ2J2 ^°m;j,-mC^,j,- M \L, 0) (£', 0| . (6) 

LL< M 

Due to the orthogonal relation of Clebsch-Gordan 
coefficients 

E r ,LM fih'M _ r f 7 \ 

° J,m;J,M-m L/ J,m;J,M-m ~ °L,L', { < ) 

m 

the on-site interaction can be reduced as the sum 
of the irreducible tensors, i.e., V — J2l = 
U^ L \(JJ)L,0)((JJ)L,0\. 

Therefore, we have proved that the engineered Hub- 
bard Hamiltonian can be written as the direct sum of 
N irreducible sub-Hamiltonians = H ( Q L) + W [L] 

= 2J X + W^ L \ The model described by each can 
be inversely mapped into a new Bloch electron model 
with a CPB, whose Hamiltonian is 

L-l 

r(L) = Yl (*i°5 a i+i + h - c -^ + ^0^0, (8) 
i=-L 

where Oj is the creation operator of new fermion and 

tj = yJJL + j + 1)(L - j). As illustrated in Fig. lc, the 
on-site interacting qubit network is reduced into a direct 
sum of N Bloch electron models with CPB. 

Spectrum- symmetry matching for nonzero U. To see 
whether the engineered Hubbard model can serves as a 
quantum data bus to coherently transfer quantum infor- 
mation, we need to study the influences of the on-site 
interaction on the SSMC. 

Let us first recall the on-site interaction free case, 
where the matrix representation of the Hamiltonian (1) 
in single-particle subspace is equivalent to that a high 
spin with angular momentum J = (N — l)/2 precessing 
in transverse magnetic field. For two-electron case, the 
above analysis shows that the original Hamiltonian can 
be reduced into N single-particle engineered models on 
Z-site chain (I = 1,3, ...,2N-1) without CPB. Obviously, 
all the eigenstates also meet the SSMC that guarantees a 
perfect QST. It is crucial for our analysis that the mirror 
symmetry is not broken even in presence of the on-site 
interaction, then the effect of nonzero U on the deviation 
of the energy levels determines the fidelity of QST via 
such system. 

Now we consider the Hamiltonian (JHJ with nonzero U . 
Since [H^ L \ R] = 0, the eigenstates can be classified into 
two sets with different parities. When U is switched on, 
a set of levels reminds unchanged while other set of levels 
deviates by a nonzero values Am — \Em(U) — Em(0)\, 
(see Fig. 2). To prove this generally, we calculate the 
action of on the eigenvectors 

L, M(~)) = e^ \L, M) = J2 d L M , M {\) \L, M>) (9) 

M' 
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FIG. 2: (color on line) (a) Numerical simulation of the energy 
levels effected by U for H {L) with L = 5, U = ~ 20. (b) 
The level shift Am of H^ L ^ for L — 20 system. It shows 
that for small M , Am is approximately uniform. Notice that 
for the optimal U = 40.5, the level shifts for small M are 
approximately equal to the half of the level difference with U = 
0. 
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This means that states \L, M(%/2)) [L — M is odd) are 
also the eigenstates of H^ L \ Indeed these corresponding 
levels are free of the on-site interaction. 

Another set of eigenvalues should be shifted by the on- 
site interaction. Imagine that if the energy deviation A m 
is not so sensitive to M, the shifts of the levels are ap- 
proximately the same. There may exist an appropriate 
U to ensure that the final levels satisfy the SSMC with 
another greatest common divisor. It will result in the 
perfect QST in the invariant subspace : {\L, M) , 

M = L, L-l,--- ,-L}. 

In order to verify our conjecture, numerical simulation 
is employed for small size systems. Exact diagonalization 
results for H^ L ' on L lattice are plotted in Fig. 2. In Fig. 
2a, the energy levels as functions of U for L = 5 system 
shows that the spectrum of consists of two sets of 
energy levels, one is independent of U, while the other is 
shifted by the repulsion. This conclusion is in agreement 
with the above analysis for arbitrary L. Numerical cal- 
culation for Am of with U = 40.5 on 21-site lattice 
is plotted in Fig. 3b. It shows that for small M, Am 
is approximately uniform. On the other hand, numerical 
calculations also indicate that in the invariant subspace 



V^ L \ the components of the state \L, — L) on the basis 
V [L ^ for small M are dominant, i.e, the effective levels of 
such state should be shifted uniformly. In other words, 
for such kind of initial state \L, —L), when U takes an ap- 
propriate value, the effective levels can satisfy the SSMC 
approximately. Thus state \L,—L) can be transferred 
into \L, L) near perfectly. 

Near-perfect swap of two electrons. Now we consider 
the spin state swapping of two electrons located on the 
two ends of the one-dimensional lattice. For the initial 
state 

|V(o)> = c\ A c^ NA |o) = I j, -j- j, —j) , (ii) 

the quantum state swapping is a mapping from \ipj(0)) 
to 



m))=C^ t Cl l \0) = -\J,J;J,J) 



(12) 



can be approximately realized in a dynamic process since 
the SSMC can be satisfied for an appropriate U as dis- 
cussed above. 

In order to confirm the above prediction, the numerical 
simulation is performed for the swapping fidelity 



F{U,t) = \ {J,J;J,J\e- lHt \J,-J;J,-J) 



(13) 



where H is the Hamiltonian (1) for the engineered Hub- 
bard model. During the time range t s [0, 10], the max- 
ima of the fidelity, F max ([/) = max{F([/, t), t ^ 10} are 
plotted in Fig. 3a as the functions of the on-site in- 
teraction strength U for the system of N = 4, 5 and 6 
sites. It shows that there indeed exist some U to get 
very high F max (/7), which seems like that the propaga- 
tion of electrons is scattering-free. One can see for the 
same N that F max have several regular peaks, and for 
each peak, the revival times T r are different. It is inter- 
esting to find that, for the first regular peak, the revival 
time T r \ ~ 7r, the second peak, T r 2 — 1.5tt, and the 
iVth peak, T rn ~ 0.57r(n + 1). The corresponding opti- 
mized U and revival time T r for each peak are listed in 
Table 1. 



Peaks U F„ 



Tr 



1 6.6 

2 11.6 

3 16.2 



0.9847 
0.9768 
0.9724 



20.6 0.9698 



3.14 ~ 7T 

4.71 ~ 3tt/2 

6.28 ~ 2tt 
7.85 ~ 5tt/3 



25.0 0.9683 9.42 ~ 3tt 



Table 1 

Table 1. The maxima of fidelity F max , the corre- 
sponding optimized U and revival time T r obtained by 
numerical simulations for the engineered Hubbard model 
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on A-site system. 

On the other hand, for different N — 2 ~ 10, more 
detailed data, such as the optimized U , F max and the 
corresponding revival time T r for the first and second 
peaks are also obtained numerically to reveal the hidden 
relationship between them. In Table 2, these numerical 
results are listed with some obvious characters for the 
cases of N — A — 10. 





1st Peaks 


2 


nd Peaks 


N 


U 


F T 

± max 1 r 


U 


max 


T 

r 


2 


2.3 


0.9999 2.72 


3.6 


0.9999 


3.50 


3 


4.9 


0.9926 3.18 


8.0 


0.9929 


4.74 


4 


6.6 


0.9847 3.14 


11.6 


0.9768 


4.71 


5 


8.6 


0.9873 3.14 


15.0 


0.9802 


4.71 


6 


10.6 


0.9906 3.14 


18.4 


0.9856 


4.71 


7 


12.6 


0.9931 3.14 


21.8 


0.9894 


4.71 


8 


14.5 


0.9948 3.14 


25.2 


0.9920 


4.71 


9 


16.5 


0.9960 3.14 


28.7 


0.9938 


4.71 


10 


18.5 


0.9968 3.14 


32.1 


0.9950 


4.71 



Table 2 

Table 2. The maximal fidelities and the corresponding 
revival time T r of the first two peaks obtained by numer- 
ical simulations for the systems with N and U . 

As for the quantum swapping scheme we study here, 
the following prominent characters can be found from the 
above analysis. Firstly, the optimized U is linearly pro- 
portional to N with slopes approximately 2 and 3.4 for 
the first and second peak respectively. Secondly, for a 
given optimized U, the larger N is, the higher the max- 
imum of the fidelity F max becomes. Consequently, the 
QST will be better. Thirdly, as the statements above, 
the revival times are tt and 37r/2 for the first and sec- 
ond peak, approximately. Our further numerical results 
shows that the above experiential law can hold for larger 
N. For different peaks, the linear relations between op- 
timized U and the sizes N are also shown in Fig. 3b. 

From the Tables we set above, it is obvious that the 
revival times obey the experience formula T rn = 0.57r(n+ 
1) approximately, where n = 1,2,..., is the order of the 
peaks of -Fmax- This formula can be understood based 
on the above analysis. If the effective levels are shifted 
by U uniformly from Am = to 1, the possible greatest 
common divisors meet the SSMC are 2/(n + 1) (n — 
1, 2, ...,). The corresponding level shift is 2[1 — l/(n+l)], 
which results in the revival period T rn . A question to be 
asked is about the relationship between the level shifts 
and the repulsion U . In Fig. 3b, the optimized U as the 
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FIG. 3: (a) Numerical simulation of the maximal -F max as 
the function of the interaction strength U for N = 4 sites 
system, (b) The optimal U as the functions of the sizes N 
obtained by the numerical simulations. The corresponding re- 
vival times, T r — tt (solid square), tt/2 (triangle), 2n (circle), 
5ir/2 (solid triangle), and 3tt (square). It indicates that the 
optimal interaction strengths U are directly proportional to N 
approximately. 



functions of the sizes N are plotted. Interestingly, they 
are simply linear functions in the range we concerned. 

Summary. In summary, we study the QST in the en- 
gineered Hubbard model with on-site interaction analyt- 
ically and numerically. It is proved that the system of 
two electrons with opposite spins in this N sites quan- 
tum network can be reduced into N one dimensional en- 
gineered single Bloch electron models with CPB rigor- 
ously. Analytic calculation and numerical results both 
show that the engineered Hubbard model can perform 
perfect quantum swapping between two distant electrons 
with opposite spins even when the certain Coulomb inter- 
action exists. The angular momentum reduction method 
could be expected to work for the QST in an engineered 
quantum spin models with more electrons. 
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